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One of the frontiers of QCD which are intensely investigated in high energy experiments is 
the high energy (small x) regime, where we expect to observe the non-linear behavior of the 
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. theory. In this regime, the growth of the parton distribution should saturate, forming a Color 

' Glass Condensate (CGC). In fact, signals of parton saturation have already been observed both in 

<^ ' ep deep inelastic scattering at HERA and in deuteron-gold collisions at RHIC. Currently, a global 



description of the existing experimental data is possible considering different phenomenological 



Qh! saturation models for the two processes within the CGC formalism. In this letter we analyze the 

1^ ' universality of these dipole cross section parameterizations and verify that they are not able to 



describe the HERA and RHIC data simultaneously. We analyze possible improvements in the 
parameterizations and propose a new parametrization for the forward dipole amplitude which 
allows us to describe quite well the small-x ep HERA data on F2 structure function as well as the 
dAu RHIC data on charged hadron spectra. It is an important signature of the universality of the 
saturation physics. 
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In the past few years much theoretical effort has been devoted towards the understanding 
of the high energy hmit of the strong interaction theory (for recent reviews see, e.g. 

In the high energy limit, perturbative Quantum Chromo dynamics (pQCD) predicts 
that the small- a; gluons in a hadron wavefunction should form a Color Glass Condensate 
(CGC), which is described by an infinite hierarchy of the coupled evolution equations for 
the correlators of Wilson lines 

mm 

. In the absence of correlations, the first equation 
in the Balitsky-JIMWLK hierarchy decouples and is then equivalent to the equation derived 
independently by Kovchegov within the dipole formalism [8]. The Color Class Condensate 
is characterized by the limitation on the maximum phase-space parton density that can be 
reached in the hadron wavefunction (parton saturation), with the transition being specified 
by a typical scale, which is energy dependent and is called saturation scale Qs [Q^ oc A^x'^]. 
Moreover, in the CGC formalism the dipole-target forward scattering amplitude J\f for a 
given impact parameter 6, which is directly related with the two-point function of Wilson 
lines, encodes all the information about the hadronic scattering, and thus about the non- 
linear and quantum effects in the hadron wave function. The function A/" can be obtained 
by solving the Balitsky-JIMWLK evolution equation in the rapidity Y = ln(l/x) [9]. Its 
main properties are: (a) for the interaction of a small dipole (r ^ 1/Qs), A/'(r) ^ r^, 
implying that this system is weakly interacting; (b) for a large dipole (r ^ 1/Qs), the 
system is strongly absorbed and therefore A/'(r) ^ 1. This property is associated to the 
large density of saturated gluons in the hadron wave function. Another remarkable feature 
of CGC formalism is that the dense, saturated system of partons to be formed in hadronic 
wave functions at high energy has universal properties, the same for all hadrons or nuclei. 

In the CGC formalism the description of the observables is directly related to the behavior 
of A/". For instance, the F2 structure function is probed in ep{A) process and is given by 

Ff(^)(x,g2) = (g2/4^2^^^)(^7MA) ^^7>(A))^ ^^^^^ Q 

/ 

= Y.j dz(fr\^>^l^^{z,r,Q^)\'' x 2 / d^hM:F{x,r,h) , (1) 

/ 

with Mj: being the fundamental representation of the forward dipole amplitude, r defining 
the relative transverse separation of the pair (dipole) and z (1 — z) the longitudinal mo- 
mentum fraction of the quark (antiquark). The photon wave functions ^ l,t are determined 
from hght cone perturbation theory (See e.g. Ref. illj]). It is useful to assume that the 



impact parameter dependence of A/V can be factorized as 7Vjr(a;, r, 6) = 7Vjr(x, r)S'(6), so 
that <7dip{x, r) = aoAfj^{x, r), with o"o being a free parameter related to the non-perturbative 
QCD physics. Similarly, the single- inclusive hadron production in hadron-hadron processes 



is described in the CGC formalism by 
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where pt and xp are the transverse momentum and the Feynman-x of the produced hadron, 
respectively. The variable Xp denotes the momentum fraction of a projectile parton and b 
is the impact parameter. Moreover, f{xp, Q'j) is the projectile parton distribution functions 
and D{z,Q'j) the parton fragmentation functions into hadrons. These quantities evolve 
according to the DGLAP [l^ evolution equations and respect the momentum sum-rule. In 
Eq. ([2]), A/jc-(fc, 6) and My^{k,b) are the fundamental and adjoint representations of the 
forward dipole amplitude in momentum space. The amplitudes N'A,j^{k, b) and Nji^riji b) 
are directly related by a Fourier transform. 

The search of signatures for the parton saturation effects has been an active subject 
of research in the last years. In particular, it has been observed that the HERA data 
at small x and low can be successfully described with the help of saturation models 
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19|]. Moreover, experimental results for the total [20||, diffractive [21 1 



and inclusive charm cross sections 
the other hand, the observed 
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231 ] present the property of geometric scaling. On 
2j] suppression of high pp hadron yields at forward rapidities 
in dAu collisions at RHIC had its behavior anticipated on the basis of CGC ideas [25|. A 
current shortcoming of these analyzes comes from the non-existence of an exact solution 
of the non-linear equation in the full kinematic range, which implies the construction of 
phenomenological models satisfying the asymptotic behavior which is under theoretical 
control. Several models for the forward dipole cross section have been used in the literature 



from Refs 
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have been able to describe the dAu RHIC data. An important aspect 



should be emphasized at this point. Although at HERA it is possible to probe values of x 



two orders of magnitude smaller than at RHIC, the saturation scales 



br these two scenarios 



are very similar due to the nuclear medium (See Fig. 1 in Ref. 27|). Consequently, one 
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can expect to be possible to cross relate these experiments in this respect and gain a clear 
understanding of the CGC in high energy experiments. There are several similarities among 
the phenomenological models proposed in Refs.[l2, 14, 1^, 
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26| . In particular, in these 



models the function M has been modeled in terms of a simple Glauber-like formula 

1 



A/'(a;, r) = 1 — exp 



(3) 



where 7 is the anomalous dimension of the target gluon distribution. The main difference 
comes from the predicted behavior for the anomalous dimension, which determines the 
transition from the non-linear to the extended geometric scaling regimes, as well as from the 
extended geometric scaling to the DGLAP regime. A detailed comparison has been presented 



in Ref. 1271. As the models from Refs. 



14 



15 



17( 1 have been exhaustively discussed in 



;he literature, in this letter we only present a brief review of the models proposed in Refs. 



12, l26(]. In the KKT model 26|| the expression for the quark dipole-target forward scattering 



amplitude is given by |26|: 

A/>(r,x) 



exp 



(4) 



where Q\ 



Cp_ 



and the anomalous dimension '~f{Y, r^) is 



, 1 + 

2 \ ^(r,r2) 



2e(F,r2) + 7C(3)c, 



(5) 



with c a free parameter ( which was fixed in 26|] to c = 4) and 



ln[l/(r2g 



sOJ 



(A/2)(F-Fo)' 



(6) 



The authors assume that the saturation scale can be expressed by Ql{Y) = A^A^/^ (^^ 
The form of the anomalous dimension is inspired by the analytical solutions to the BFKL 
equation 28|]. Namely, in the limit r — with Y fixed we recover the anomalous dimension 
in the double logarithmic approximation 7 ^ 1 — ^l/{2^). In another limit of large Y with r 
fixed, Eq. (JSj) reduces to the expression of the anomalous dimension near the saddle point in 



the leading logarithmic approximation 7 



~ i _|_ 

2 ^ 14cC(3) 



. Therefore Eq 



jnj mimicks the onset 



of the geometric scaling region 17|, 129(1 . In the calculations of Ref. 26|] it is assumed that a 
characteristic value of r is r ^ 1/(2 k^) where is the transverse momentum of the valence 
quark and 7 was approximated by 7(F, r^) ^ 7(y, 1/(4 In the above expressions the 



parameters A = 0.6 GeV and A = 0.3 are fixed by DIS data [l^. Moreover, the authors 
assume Iq = 0.6. The initial saturation scale used in is defined by QIq = Q1{Yq) 
with Yq being the lowest value of rapidity at which the low-x quantum evolution effects 



are essential. As demonstrated in Ref. 



26| this parameterization is able to describe the 



dAu RHIC data when the forward dipole cross section is convoluted with the respective 
fragmentation function and the parton distributions for the deuteron. On the other hand, 
in Ref. [l^ another phenomenological saturation model has been proposed in order to 
describe the dAu RHIC data (hereafter denoted DHJ model). The basic modification with 
respect to the KKT model is the parameterization for the anomalous dimension which is 
now given by 

7(r,r2) =7. + A7(r,r2) (7) 



\l0g j;T0t\ 



where 

Arv(V = n - 

XY + \hgj;^\+dVY 

with Qt = Qs{Y) a typical hard scale in the process, A = 0.3 and d = 1.2. Moreover, 7^ = 
0.63 is the anomalous dimension for BFKL evolution with saturation boundary condition. 
Similarly to the KKT model this model is able to describe the dAu RHIC data. 

As already discussed in Ref. |27|] , based on the universality of the hadronic wave function 
predicted by the CGC formalism, we might expect that the KKT and DHJ parameteriza- 
tions would also describe the HERA data on proton structure functions in the kinematical 
region where the saturation effects should be present (small x and low Q^). However, this 
expectation fails when the KKT model is applied, as verified in [27]. Here we extend the 
analysis to the DHJ model without any modification of the parameters fitted at RHIC, 
only assuming A = 1 and adjusting the non-perturbative parameter ctq, which defines the 
normalization, in order to describe the F2 experimental data at = 10 GeV^. In Figs. [1] 
and [2] we present the predictions of the DHJ model for the proton structure function and 
compare with the ZEUS data 301]. We can see that this parameterization fails for both small 
and large values of Q^. Consequently, the current parameterizations of the forward dipole 
cross section which are constrained at RHIC are not able to describe the HERA data. An 
open question is if minimal modifications in these parameterizations allow to describe both 
sets of data. Following Ref. |3l| we consider a modification of the KKT model assuming 
that the saturation momentum scale is given as in the GBW model, Yq = 4.6, c = 0.2 and 
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FIG. 1: The proton structure function at different values of the photon virtuaUties. Data from 
ZEUS. 

that the typical scale in the computation of ^(Y, r^) is the photon virtuality. Its predictions 
(KKTm lines) are presented in Figs. [T]and[2l It is observed that these modifications imply 
a quite good description of the HERA data. Similarly, as the evolution of the F2 data is 
not well described by the DHJ model it is possible to improve this model by the modification 
of the anomalous dimension. Here we propose to modify the DHJ model assuming now that 
Qt = Qo = 1-0 GeV, i. e. that the typical scale is energy independent. It is important to 
emphasize that this modification preserves the main properties of the anomalous dimension 
proposed in [12]. Basically, we still have that the anomalous dimension increases logarith- 
mically with pt from 7 = 7^ to its asymptotic value 7^1, while decreasing with Y as 
A7 ^ 1/Y at very large rapidity. As shown in Figs. [T] and [21 with this modification our 
predictions (GKMN lines) agree with the experimental data. 

The question which follows is whether the RHIC data are still well reproduced after these 
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FIG. 2: The proton structure function at different values of the photon virtualities. Data from 
ZEUS. 



modications. Following Ref. 12|] we have calculated the single inclusive hadron production 



cross section in dAu collisions at different rapidities. We have used the CTEQ5L quar k and 



gluon distributions [32] and the LO KKP quark-hadron fragmentation functions 33|]. Our 



results are presented in Fig. [3] and compared with the BRAHMS data 2J]. The KKTm and 



GKMN predictions are represented by long-dashed and solid curves respectively. As in Ref. 



12l | we need a i^-factor in our calculations, since it has been performed at leading order in a^. 
Although the normalization should be modifed by these corrections we expect that the shape 
of the momentum distributions should not change. Our values of K have been determined 
so as to reproduce the data at pt = 1.0 GeV and they depend on the parameterization 
adopted. For KKTm we find a larger value of K than for GKMN. Moreover, while the KKTm 
parameterization fails to describe the full set of data, the GKMN one is able to reproduce 
the data quite well even at very small values of px- Consequently, the GKMN model is able 



7 



> 

0) 

C5 



10^ 

10" 



Q. 

CM 

■o 
>« 
■o 



r 10 r 



10" 

10"' 



TJ 10" 

10"' 




^GKMN ^'^ 

=5.0 

KKTm 



10^ 
10" 

-t 10"' t 
10" 
10"' 



10"- 



I . I i I . I i I . I i I . I i I 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



> 
0) 



10 

10° 
10"^ 



Q. 

CM 

■o 
>« 
■o 



r 10 



10" 
10"' 



TJ 10"- 

10"' 



I ' I ' I ' I ' I ' I ' I ' I ' 
y = 2.2 -: 




I . I i I . I i I . I i I . I i I 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

P,[GeV] 



10 



M ' I ' I ' I ' I ' I ' I ' I ' I ' 

y = 1 - 



I . I i I . I i I . I i I . I i I 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



10 
10° 
10"^ 
10"^ 
10"^ 
10"^ 
10"^ 
10"^ 



' I ' I ' I ' I ' 

y = 



3.2 J 








0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

P,[GeV] 



FIG. 3: Comparison of theory and BRAHMS data for minimum-bias dAu collisions at RHIC 
energy. 

to describe the ep HERA and dAu RHIC data in terms of an unique parameterization for 
the dipole scattering amplitude, which is based on the saturation physics. 

Before presenting a summary of our main results, let us briefly discuss the basic properties 
of the resulting GKMN model (a more detailed analysis will be presented elsewhere). In Fig. 
4 (a) we present the forward dipole cross section as a function of the scaling variable rQs 
for distinct parameterizations. As it can be seen the DHJ, KKTm and GKMN models have 
a similar behavior. The difference among the models can be demonstrated studying the 
behaviour of the effective anomalous dimension, defined by 7e// = '''di^Xq^i/l-^'^ (See similar 
analyzes in Ref. 3l|)- In Fig. 4 (b) is shown 7e// as a function of the virtuality Q^, using the 
average dipole size as r = 2/Q. While the GBW model presents a fast convergence to the 
DGLAP anomalous dimension at large Q^, the IIM parameterization has a mild growth with 
virtuality, converging to 7 ~ 0.85 at large Q^. The KKTm and IIM parameterizations are 
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similar at large Q^, but differ at small virtualities, with the KKTm one predicting a smaller 
value. On the other hand, the predictions of the DHJ and GKMN parameterizations are 
similar at small and differ at large virtualities. In particular, we have a strict difference 
between these models in the intermediate range of virtualities, which can explain why the 
DHJ model does not describe the evolution of the F2 structure function. 

As a summary, in this letter we have analyzed current parameterizations for the dipole 
scattering amplitude which are able to describe separately the ep HERA and dAu RHIC data. 
We have shown that an unified description using these parameterizations is not possible. 
We have proposed a modification in the DHJ parameterization for the dipole scattering 
amplitude, based on saturation physics, which allows to describe simultaneously the ep 
HERA and dAu RHIC data. This result has been obtained adjusting the normalizations of 
the dipole cross section and single inclusive hadron cross section and assuming an energy 
independent typical scale, keeping all other original parameters. A global least fit of data 
would change slightly the values of our parameters. This would be a fine tunning which 
is beyond the scope of this work. We rather prefer to keep the level of fitting accuracy 



of and emphasize the strategy to reconcile two different sets of data. Apart from 
this fine tunning, a more detailed theoretical study of the proposed anomalous dimension 
is necessary. We postpone these improvements for a future publication. Finally, our results 
demonstrate that an unified description of the experimental data which probes the high 
energy regime of QCD is possible. This is an important signature of the universality of the 
saturation physics. 
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